Dissipative dynamics of quantum discord under quantum chaotic environment 
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We investigate the dissipative dynamics of quantum discord in a decoherence model with two 
initially entangled qubits in addition to a quantum kicked top. The two qubits are uncoupled 
during the period of our study and one of them interacts with the quantum kicked top. We find 
that the long time behavior of quantum discord could be well described by the fidelity decay of the 
quantum kicked top; for short time behavior, however, the phase of the amplitude of the fidelity 
decay is necessary to provide more specific information about the system. We have made comparison 
between the quantum kicked top and multi-mode oscillator system in describing environment, and 
also compared the dynamics of the entanglement with that of quantum discord. 
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It was a surprise that some quantum computing ap- 
proaches could be carried out efficiently in the absence 
of entanglement [l[ . This implies that entanglement does 
not involve complete quantum correlation and thereby 
a new terminology 'quantum discord' (QD) has drawn 
more and more attention [2|-|5( recently and has been 
considered to replace entanglement as a vital quantum 
resource for quantum information processing. 

Since no quantum system is isolated from the envi- 
ronment, it is of fundamental and practical importance 
to study the dissipative dynamics of the QD. Modeling 
the environment by a multi-mode oscillator system [6[ 
had been widely adopted, but recent research pointed 
out that the environment could also be simulated by 
the quantum chaos model Q, such as a quantum kicked 
top (QKT) Q. It was shown that the quantum chaos 
model could simplify the treatment about environment 
compared with the multi-mode oscillator system model, 
and also efficiently simulate both Markovian and non- 
Markovian environments 0] . 

Our focus in the present work is on the dynamics of 
the QD under a quantum chaos model, in which the 
decoherence is from a QKT and the decoherence 

could be reflected by fidelity decay (FD). We consider 
two qubits with one of them coupled to a QKT. The two 
qubits are initially entangled, but with no coupling dur- 
ing the period of our investigation. Our purpose is to 
investigate the physics behind the dynamics of the QD 
under this QKT-induced decoherence. To this end, we 
will check the dynamics of the QD with respect to the 
chaotic and regular regimes of the QKT. We find that 
FD is suitable for describing the long-term dynamics of 
the QD, which shows difference between chaotic and reg- 
ular QKTs. Moreover, from the phase of the amplitude 
of FD (PAFD), we may discover more differences in the 
environment-induced behavior between chaotic and reg- 
ular QKTs. Furthermore, we will try to explore the simi- 
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larity, from the dynamics of the QD, between the chaotic 
(regular) QKTs and Markovian (non-Markovian) multi- 
mode oscillator environment. In addition, it has been 
shown that the relative entropy entanglement (REE) is 
closely related to the QD [llj. As a result, a comparison 
of the QD with the REE would make us understand the 
dissipative dynamics from another angle. 

Let us first briefly review the main points of QD. 
The QD is defined to quantify the difference of mu- 
tual information (MI) between the classical and quantum 
forms. In classical information theory, the correlation 
between two random variables A and B is described by 
MI using following two expressions: I(A, B) = H(A) + 
H(B) - H(A, B) and J(A, B) = H(A) - H(A\ B), where 
H(A) and H(B) are the Shannon entropies, H(A,B) is 
the joint Shannon entropy, and H(A\B) is the condi- 
tional entropy introduced for quantifying the ignorance 
about the value of A given a known B. The two ex- 
pressions above are equivalent in the classical case using 
Bayes' rule, and their quantum version is the replacement 
of the Shannon entropy by von Neumann entropy, i.e., 
I{pa,Pb) = S(p A ) + S(p B ) - S{pa,b), and J(pa,Pb) = 
S{pa) — S(pa \ Pb), where quantum conditional entropy 
S(pa\pb) depends on the choice of measurement. If 
we perform a set of projectors {IF^} on the subsys- 
tem B, then the post-measurement states of the sub- 
system A with outcomes j is Pa\h b = PjTr[Hf pAB^lf] 
with pj = Tr[I <E) Pab], and the quantum conditional 
entropy is defined as S(pa\pb) = J2PjS(.PA\n B )- So 
different choices of {n^} result in different S(pa\pb), 
and the QD is defined as the minimum difference of 
the quantum MI between I(pa,Pb) and J(pa, Pb), i-e., 
Q = min[I(pA,PB) - J{pa,Pb)}, where J{pa,Pb) is also 
called classical correlation (CC) |12j . 

Consider the QKT with following Hamiltonian [Io| 

oo 

He = t Jx + 2j j * S ^-^n a) 

n— — oo 

where J x and J z are angular momentum components 
with J = J x i + J y j + J z k, T is the time period, and 
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FIG. 1: (color online) Time evolution of the FD for chaotic 
(blue solid line) and regular (red dot-dashed line) QKTs in 
units of T, where the parameters are J = 100, v = n/2 and 
e = 0.001. We set r\ = 20 and 0.1, respectively, for chaotic 
and regular QKTs. The inset shows the time evolution of 
the PAFD a for chaotic (blue solid line) and regular (red 
dot-dashed line) QKTs in short time scale, where, since the 
variation in the regular case is much larger than that in the 
chaotic case, we enlarge a by 40 times in the chaotic case, for 
a comparison with the regular case. 



v and rj represent the precessional angle and the kicking 
strength, respectively. We assume that both the dynam- 
ics of the qubits and dissipations are much slower than 
dephasing. So we may only focus on the dephasing effect. 
We consider one of the qubits (A) is free and the other 
(B) interacts with QKT via the Hamiltonian 



H = eafV E + H E 



(2) 



where e is the strength of the interaction between the 
qubit and the QKT, and af is the Pauli operator of the 
qubit B. V E is an operator of the QKT, which, for conve- 
nience, could be chosen as J x in the following calculation. 
The corresponding Floquet operator can be constructed 
as 



(3) 



In the classical limit J — » +oo, the chaoticity degree is 
determined by the kicking strength 77 and the precessional 
angle v. For example, for v — n/2, the classical kicked 
top is dominated by a regular motion for 77 < 2.5, whereas 
the motion is prevailed by chaos for 77 > 3 [13]. The 
regime with 2.5 < 77 < 3 is for mixture of chaos and 
regularity. For our purpose in this work, we will only 
consider the regular and chaotic regimes, instead of the 
mixture regime. 

Consider following initial state of the system, 



Pabe = 



3=x,y,z 



>4>\, (4) 



where Cj are real numbers, and (&f) is the usual 
Pauli matrix acting on the subsystem A (B) 2]. The 



initial QKT state \6,, 
as [iH 



is a spin coherent state defined 



\6,4>) = exp[—i8(J x sin <j) — J y cos </>)] | J, J) , (5) 

where | J, J) is the eigenstate of J 2 and J z with eigenval- 
ues J (J + 1) and J, respectively. The time evolution of 
the density matrix after n kicks can be described as 



Pabe = u n PABE {u^Y 



(6) 



with the corresponding reduced density matrix of the two 



qubits, in the basis {|1) = |0) A |0) B , |2) 
|3) = |1) A |0) B , |4) = |1)J1) B }, given by 
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(0,4>\ exp[i(eV E 

is the amplitude of 

which is generally complex, i.e., 
2 . 



where / = (#'(t)|#(i)) 

H E )t]exp[-i(-eV E + H E )t] 
FD of the QKT [F 

/ = VFe la , where F — is the FD, and the PAFD a 
reflects the relative phase between the QKT states \^(t)) 
and \^/'(t)). For a fully chaotic quantum system real- 
ized through adjusting the parameter 77, FD behaves with 
Gaussian or exponential decay depending on the values 
of the perturbation strength [15J . Change of the param- 
eter 77 moves the quantum system to the regular regime 
of the QKT, and then yields only the Gaussian decay 
and revivals of the FD (15l [l6| . The corresponding re- 
vival time is r — k/e with k the constant. Fig. 1 shows 
the FD for chaotic and regular QKTs. In our calcula- 
tion here and following, the time unit is the period T, 
and the initial state \6, 4>) is a randomly chosen spin co- 
herent state. As shown in Fig. 1, the fluctuation in the 
chaotic case could be understood as being from the fi- 
nite numbers of the eigenstates of the Flo quet operator 
U involved in the initial QKT state \9,4>) [131, an d the 
revivals in the regular case is related to the underlying 
periodic classical motion [l6j]. Both the fluctuation and 
the revival could be regarded as the memory effects. The 
inset in Fig. 1 shows the regular and chaotic behaviors 
reflected by PAFD: Their variations are regular and ir- 
regular respectively, and the amplitude of the variation 
in the regular case is much larger than that in the chaotic 
case. 

The reduced density matrix of the two qubits in Eq. 
(7) is an X state, whose QD and CC could be calculated 
analytically, using the methods in [l8j], as 



Q{p'a,b) = 2 + log 2 Xi - C[p' AtB ), (8) 



and 
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C(p^ B ) = l-min[-^ 
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log. 



1 + (-l) k 6 2 , 



(9) 



where Ai = 

c z - \(cx - Cy\, 
A 4 = [1 - c z - 



C Z + \(C X - Cy)\ 

F]/4, A 3 = [1 - c z 
{c x +c y )\ VF]/4, 



A 2 = [1 + 
,+ Cy )|/F]/4, 
IcJ, and 



72 = 

[2(4 + eg) + 2 | C 2 - eg | (|cos2a| + |sin 2a\)]F/2. For 
a real or imaginary amplitude of the FD, 2 could be 
reduced to 6> 2 = yF max[|Ca.| ,\c y \], which is similar to 
the result of the phase damping channel in Q, describ- 
ing the dephasing induced by the Markovian multi-mode 
oscillator environment. 

To be more clarified, we compare the QD with the REE 
for entanglement . Direct calculation yields the REE 
in our case, 



E = l + log 2 p + (1 - /?') log 2 (l - /?'), 



(10) 



where /3' = 
max[Ai, A 2 , A 3 , A4 



max [0.5, A max ], with A max = 
We plot the evolution of QD 
and REE, along with CC and FD, in Figs. 2 and 3 for 
a comparison. In the chaotic regime, both the QD and 
REE are only determined by the FD, but in the regular 
regime, the REE still only depends on the FD, whereas 
QD behaves complicated: Its long time behavior is fully 
determined by FD; the short time behavior is more 
relevant to the PAFD a. a. 

To further check the Markovian and non-Markovian 
dephasing effects of the QKT, we consider below the 
correspondence associated with the chaotic and regular 
QKTs, respectively. For the chaotic QKT, the PAFD a 
is nearly zero during the short time evolution of the FD. 
So the effect from the chaotic QKT is similar to the phase 
damping channel, a typical Markovian environment. Like 
the model of the phase damping channel, there are three 
kinds of CC dynamics in our case here, depending on the 
parameters Cj, i.e., the constant CC, CC with a sudden 
change, and CC with monotonic decay 0. Our interest 
is in the case with a sudden change, where the QD can 
be maintained for some time [20j . To show similarities 
and differences between our decoherence model and the 
phase damping channel, the dynamics of QD, CC and 
others, whose short time behaviors are very similar to 
those under the phase damping channel [2(|, are plotted 
in Fig. 2. But for long time behavior, the situation is 
different: For example, there are fluctuations in the evo- 
lution of the QD for a finite angular momentum J, which 
does not happen in the phase damping channel. In ad- 
dition, we could find from the Fig. 2 (a) the consistency 
in the short time evolution between the CC and the FD, 
and also between the REE and the FD. Furthermore, in 
comparison to the sudden death in the time evolution of 
REE, there is only instantaneous disappearance of the 
QD at some time points as shown in Fig. 2 (b), where 
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FIG. 2: (color online) (a) Short time evolution of the QD (blue 
dashed line), CC (red dash dotted line), REE (green solid 
line), and FD (black dotted line) in units of T with respect 
to the chaotic QKT, where the parameters of the FD are the 
same as in Fig. 1 and the initial state is with the parameters 
c x =0.95, c y =-0.85, and c z =0.85. (b) Long time evolution 
of the QD, CC, REE, and FD under the same parameters in 
(a). The inset is a magnification for the time period 2900 < 
t/T < 3000. 



the small fluctuations of the QD are originated from the 
memory effects of the QKT. Those fluctuations disap- 
pear only in the classical limit, i.e., J — > +00 in our case, 
which corresponds strictly to the Markovian regime. 

The memory effects of the QKT could be more evident 
in the treatment of the regular QKT, which resembles the 
dynamics in the non-Markovian environment [f|. From 
short time behaviors, we may find that before the CC 
becomes a constant, there exist oscillations in both the 
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FIG. 3: (color online) (a) Short time evolution of the QD (blue 
dashed line), CC (red dash dotted line), REE (green solid 
line), and the compressed quantity |cos2a| + |sin2a| (black 
dotted line) relevant to the regular QKT in units of T, where 
other parameters are the same as in Fig. 2 except r\ =0.1. 
(b) Long time evolution of the QD (red solid line), CC (blue 
dashed line), REE (green solid line), and the quantity FD 
(black dotted line) under the same parameters in (a). The 
inset is a magnification for the time period 2900 < t/T < 
3300. 



CC and the QD, as shown in Fig. 3 (a). Although these 
oscillations are very similar to those in a non-Markovian 
environment ||, their origins are different. The oscilla- 
tion in the non-Markovian environment is related to the 
dephasing (2lj . whereas in our case it comes from the 
PAFD. In contrast, there is no oscillation in the short 
time behavior of the REE (see Fig. 3(a)), which could be 
considered as an important difference between entangle- 



ment and QD. Another difference between entanglement 
and QD could be found in the long time evolution (See 
Fig. 3(b)), where REE experiences sudden death and 
birth, but not for QD. 

Therefore, only the chaotic QKT in the classical limit 
describes the Markovian environment. The rest of the 
QKT, such as the finite J of the chaotic QKT as well 
as the regular QKT, corresponds to the non-Markovian 
environment. In this sense, the results in [20I ] are only 
relevant to the chaotic QKT in classical limit. It would 
be of great interests to revisit the problems in [20j by 
regular QKT or finite J of chaotic QKT. Moreover, en- 
tanglement is usually described by concurrence. If con- 
currence is employed in our model, however, the main 
results obtained above would be remained, because both 
REE and concurrence involve no PAFD. So entanglement 
could be distinguished from QD by PAFD regarding the 
regular regime of the environment. 

Our model could be straightforwardly extended to 
some more complicated cases, such as two qubits inter- 
acting with a QKT commonly or with two QKTs inde- 
pendently, particularly from the initial state Eq. (4). 
As a result, our model would be useful to study either 
two locally located quibts experiencing the same envi- 
ronment or two distant qubits suffering different deco- 
herence. Moreover, our model involves only dephasing 
induced by the QKT. It might be more interesting to 
take into account a model involving more complicated 
dissipation in the future. 

In summary, we have studied the dissipative dynamics 
of the QD for two qubits initially entangled, but only one 
of them coupled by a QKT. We have checked the dynam- 
ics of the QD in both the chaotic and regular regimes of 
the QKT, in comparison to entanglement. A compar- 
ison of the dynamics of the QD under QKT with the 
multi-mode oscillator model has also been made, which 
would help to go beyond multi-mode oscillator approxi- 
mation and to understand dissipative dynamics from an- 
other viewpoint. Since QD could be used to understand 
the efficiency of quantum gate operations in the absence 
of entanglement and quantum systems might suffer from 
quantum chaotic environment, our study should be help- 
ful for not only further understanding quantum chaos, 
but also suppressing decoherence in quantum informa- 
tion processing. 
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